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Abstract 

We introduce a model for charge and heat transport based on the Landauer-Biittiker 
scattering approach. The system consists of a chain of N quantum dots, each of them 
being coupled to a particle reservoir. Additionally, the left and right ends of the chain 
are coupled to two particle reservoirs. All these reservoirs are independent and can be 
described by any of the standard physical distributions: Maxwell-Boltzmann, Fermi- 
Dirac and Bose-Einstein. In the linear response regime, and under some assumptions, 
we first describe the general transport properties of the system. Then we impose the 
self-consistency condition, i.e. we fix the boundary values (Tl, /xl) and (Tr, /xr), and 
adjust the parameters {Ti, jii), for i = 1, . . . , A^, so that the net average electric and 
heat currents into all the intermediate reservoirs vanish. This condition leads to ex- 
pressions for the temperature and chemical potential profiles along the system, which 
turn out to be independent of the distribution describing the reservoirs. We also de- 
termine the average electric and heat currents flowing through the system and present 
some numerical results, using random matrix theory, showing that these currents are 
typically governed by Ohm and Fourier laws. 

Mathematics Subject Classification 80A20, 81Q50, 81U20, 82C70 

Keywords Quantum transport; Quantum dots; Landauer-Biittiker scattering approach; 
Onsager relations; Entropy production; Random matrix theory; Ohm and Fourier laws 

1 Introduction 

The study of transport properties of systems out of equilibrium is a fascinating subject in 
theoretical physics. In particular, various models have been developed to find out what 
are the underlying microscopic mechanisms giving rise to macroscopic behaviours such as 
Ohm and Fourier laws [1-3]. 
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One interesting class of models that has been introduced with this purpose is related to 
the Lorentz gas [4-12]. In particular, a classical model for particle and energy transport has 
recently been investigated by Eckmann and Young [12]. In this EY-model, the system is 
a linear chain of chaotic cells, each one containing a fixed freely rotating disc at its centre 
(see Figure 1). Both ends of the chain are coupled to classical particle reservoirs assumed 
to be in thermal equilibrium at temperatures Tl and Tr, respectively. Particles are injected 
into the system from these reservoirs at characteristic rates 7l and 7r, respectively. Once 
in the system, the non-interacting particles move freely and collide elastically with the 
boundaries of the cells as well as with the discs. The particles can also leave the system 
through these reservoirs. 

Note that in this model the discs play the role of energy tanks, allowing the redistribu- 
tion of energy among the particles in a given cell and thus permitting the system to reach at 
large times a stationary state satisfying local thermal equilibrium (LTE). As a consequence, 
this model admits a well-defined notion of local temperature. 

One of the main results states that, under some assumptions, the average stationary 
particle and energy currents through the system are governed by Pick and Fourier laws. 

In this paper, we will always assume that the considered system has reached its (unique) 
stationary state and discuss only its average transport properties. In particular, we shall 
only consider average stationary currents and consequently simply call them currents. 




(1) (2) (N) 

Figure 1 : The classical EY-model composed of N cells. 

In the present work, we shall construct a similar type of model using the Landauer- 
Buttiker scattering approach [13-17]. This will permit, in particular, to establish an effec- 
tive quantum version of the EY-model. 

Basically, in the scattering approach to charge and heat transport, one expresses the 
(average) electric and heat currents through a conductor in terms of the scattering data of 
this conductor. Therefore, starting from the EY-model, we shall make the following crucial 
modifications: 

(1) The main modification concerns the scattering version of the disc. There are two 
features of the EY-model that should be pointed out: (i) The exact position of the 
disc within a given cell is not important, so one may consider the cell represented in 
Figure 2 (left), (ii) Once the system has reached a stationary state satisfying LTE, 
the boxes containing the discs are at equilibrium and consequently the net (average) 
particle and energy currents into these boxes vanish [7-9]. 
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Figure 2: Left: A possible EY-cell. Right: The corresponding scattering cell. 

Therefore, in the scattering approach, one can obtain an effective description of the 
discs as follows^: We replace the boxes containing the discs by independent particle 
reservoirs with characteristic temperatures and chemical potentials such that the net 
(average) electric and heat currents into these particle reservoirs vanish (Figures 2 
and 3). In this manner, the particles in a given cell will be in equilibrium with their 
respective reservoir, i.e. the system will be in local thermal equilibrium. 

Although the boxes containing the discs have finite energy, while the reservoirs have 
infinite energy, we believe that the transport properties of the system in the stationary 
state (assuming there is one) will be similar in both cases. 

Usually, the condition of zero net (average) particle and heat currents, between the 
system and a reservoir, is referred to as the self-consistency condition [18-23]. If 
one is only interested in the electric current, the self-consistent reservoirs correspond 
to the voltage probes used in mesoscopic physics, originally introduced to model the 
inelastic scattering occurring in a conductor [24-32]. 

A crucial distinction appears here: In voltage probes only the electric currents are set 
to zero while in this paper the self-consistent reservoirs require to set both the electric 
and heat currents to zero. To our knowledge, this fundamental difference has never 
been discussed in the formalism of Landauer-Biittiker and consequently makes one 
of the novel aspects of the present work. In particular, our results generalise works 
by, for example, Buttiker [24-26], D' Amato and Pastawski [27], Ando [30] and very 
recently Roy and Dhar [23]. 

If the couplings between the system and the self-consistent reservoirs are sufficiently 
small, then one may interpret the self -consistent reservoirs as ideal potentiometers 
and thermometers [33, 34]. 

(2) The transport properties of the cells will be given in terms of N scattering matrices 
S^^\ . . . , 5"^^^. However, we will work mostly in terms of the scattering matrix S 
associated to the global multi-terminal system, which is obtained by composing the 
local scattering matrices together. 



This way of modelling the discs in the scattering approach was suggested to me by M. Buttiker. 



4 



Philippe A. Jacquet 



(3) In order to have some generality, we will describe the particle reservoirs in terms 
of the three standard physical distributions: the classical Maxwell-Boltzmann dis- 
tribution and the quantum Fermi-Dirac and Bose-Einstein distributions. We assume 
that the particles may carry some charge e, so that we may speak of electric currents 
instead of particle currents. Nevertheless, we shall assume, as in the EY-model, that 
the particles do not interact with each other. 




Figure 3: The system made of N quantum dots with self-consistent reservoirs. 

For obvious reasons, we now call the cells quantum dots. In the EY-model, the outer 
boundaries of the cells were chosen such that the system was chaotic. We shall at some 
point impose this property by assuming that the quantum dots are classically chaotic. 

The organisation of this paper is as follows. In Section 2, we establish the framework of 
multi-terminal systems in mesoscopic physics, present our model and outline the strategy 
adopted to obtain the desired chain of quantum dots. Starting with Section 3, we will work 
in the linear response regime and under the following main assumption: the scattering 
matrix S associated to the system does not depend on the energy. As we shall see, this 
assumption, which is a good approximation in certain limiting cases, will lead to some 
interesting consequences. 

More precisely, the main results of this paper are divided into three parts. In Section 3, 
we present the general transport properties of the system, and discuss in particular the 
Onsager relations and the entropy production. In Section 4, we solve the self-consistency 
condition, i.e. we fix the boundary values (Tl, /il) and (Tr, /xr), and adjust the parameters 
{Ti, Hi), for i = 1, . . . , N, so that the net average electric and heat currents into all the 
intermediate reservoirs vanish. This condition leads to expressions for the temperature and 
chemical potential profiles along the system, which turn out to be independent of the nature 
of the particles, i.e. independent of the distribution describing the reservoirs. From these 
profiles, we work out the average electric and heat currents flowing through the system. 
Let us emphasise that all these results are obtained rigorously and actually hold for any 
geometry (see Figure 4). These results constitute the main analytical results of this paper. 

In Section 5, we restrict our attention to the linear geometry represented in Figure 3. In 
analogy with the EY-model, we assume that the quantum dots are classically chaotic and 
consequently use random matrix theory (RMT) to describe the transport properties of the 
quantum dots. While we cannot prove the validity of Ohm and Fourier laws in this situa- 
tion, some approximate derivations are feasible in some particular cases (see Remark 5.5). 
We thus turn to numerical simulations and found that Ohm and Fourier laws typically hold 
(in the sense of RMT) in a chain of quantum dots with self-consistent reservoirs. 
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2 The Model 

We consider a chain of N connected quantum dots, where each dot is coupled to a particle 
reservoir at temperature Tj and chemical potential fXi, with i — 1, . . . , N. Additionally, the 

left and right ends of the chain are coupled to particle reservoirs with parameters (Tlj/xl) 
and (Tr,/xr), respectively (Figure 3). All these reservoirs, which we also call terminals, are 
independent and inject particles into the system according to some distribution function. 
We assume that they can also absorb particles without changing their state. In this paper, 
we shall consider the following three cases (i e {L, R, 1, . . . , N}): 

fr{E) ^ r''{E;T„f,,) = exp (-^^) , (2.1) 

(2.2) 
(2.3) 

Here, E E [0, oo) is the energy, /cb is Boltzmann's constant and Tj > is the temperature 
of the i-th terminal. For the chemical potentials there is one subtlety: for the Maxwell- 
Boltzmann and Fermi-Dirac functions, one has e (— oo, oo), while for the Bose-Einstein 
function, one has G (— oo, 0). Let us emphasise that, for each i, the parameters Tj and 
/li are independent of each other. 

We consider that the transport properties of the k-th quantum dot are described by a 
scattering matrix S'^'^^ (E) at energy E > 0: 

S^'\E) = (s^l^iE)) , ke{l,...,N}, (2.4) 

where i,j denote the three possible entrances of the k-th dot and the indices m, n denote 
their corresponding channels. In Section 5, we shall use these scattering matrices to build 
the scattering matrix S associated to the global system made of N quantum dots. For the 
present, let us assume that we are given the scattering matrix of the global system: 

S{E) = {Sij.mn{E)) , (2.5) 

where i,j e {L, R, 1, . . . , A^} and for each couple the indices m E {1, . . . , Mj} 
and n E {1, . . . , Mj} number the channels in terminal i and j, respectively. Therefore, 
the matrix element Sij-mniE) is the probability amplitude that a particle with energy E 
incident in channel n in terminal j is transmitted into channel m in terminal i. 

Being given S{E), one can define the total transmission probability tij{E) that a parti- 
cle with energy E goes from terminal j to terminal i [35]: 

Mi Mj 

tij{E)^J2J2\^^r,mn{E)\'. (2.6) 

m=l n=l 

Throughout this paper, we always assume that the following holds. 



fr{E)^r^{E;T„t,,)^ 



fFD 



exp 



E - Hi 
knTi 



+ 1 



fr{E)^r^{E;T,,i,,) 



exp 



E — jjLi 

~kWI\ 



6 



Philippe A. Jacquet 



Assumption Al. For each E > 0, the complex matrix S{E) is unitary. 

On the other hand, we do not assume the scattering matrix S{E) to be symmetric. 
Using the unitarity of S{E) one easily obtains the following results: 

J2U,{E) = M, , Wj and ^ t,, (E) = M, , . (2.7) 

i 3 

Remark 2.1. The sums appearing in (2.7) are over the set {L, R, 1, . . . , N} of all terminals. 
Unless stated, in what follows every sum will be understood over this set. 

2.1 The Currents 

In this paper, we shall always assume that the following holds. 

Assumption A2. ( i) The particles do not interact with each other 

(ii) The system admits a unique stationary state (out of equilibrium). 

Under these assumptions, one can derive for any of the distributions (2.1)-(2.3) the 
expressions for the (average) electric currents in a multi-terminal system [17, 29, 36, 37]. 
One finds 

= T: E / [t3iiE)fi{E) - Uj{E)fj{E)] dE , (2.8) 

where e > is the charge carried by the particles and h is Planck's constant. 

Remark 2.2. To our knowledge, the expression (2.8) for the average electric current in a 
multi-terminal system was first introduced, through formal arguments, by Biittiker in [17] 
and was very recently given a rigorous and general derivation (together with (2.9)) in a 
work by Aschbacher et al. [37] using a C* algebraic approach. 

The expression (2.8) can be interpreted as follows: tji{E) fi{E) is the average num- 
ber of particles with energy E that are transmitted from terminal i to terminal j, and 
tij{E)fj(E) is the same but from terminal j to terminal i. Therefore, Ii is the net aver- 
age electric current at terminal i, counted positively from the terminal to the system. 

Remark 2.3. To describe the transport of neutral particles (such as neutrons or phonons), it 
suffices to set e = 1 in (2.8) and one will obtain the particle current. 

It turns out that most relations that we shall encounter will have the same form in the 
three considered cases (an example being given by the expression (2.8)). Therefore, from 
now on, all relations not wearing the superscript MB, FD or BE will be assumed to hold in 
all three cases. 

We shall also obtain some expressions which do not depend at all on the distribution 
function describing the reservoirs. In order to emphasise the fact that such relations do not 
depend on the nature of the particles, we shall say that they are universal. We hope the 
reader will not be confused with the notion of universality used to refer to properties that 
are model independent. 
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Under the Assumption A2, one can also derive the expression for the (average) heat 
current at terminal i [37, 38]. Basically, the idea is to write the particle current, which is 
given by (2.8) with e = 1, as well as the energy current, which is also given by (2.8) with 
e — 1 but with an extra E in the integrand, and then to invoke the first law of thermody- 
namics 5Qi — dEi — jiidNi. One finds 

J^^h^j^ ME)ME) - kj{E)fj{E)] {E - tx,)dE . (2.9) 

Remark 2.4. As we shall see in the linear response analysis, if one considers Sjjii/e and 
5Ti/T as the thermodynamic forces, then 7j and Jj turn out to be the right currents to 
obtain the Onsager relations. 

Although the expressions (2.8)-(2.9) for the electric and heat currents have a direct 
physical interpretation in terms of tij{E), we prefer to work in terms of the quantities 

r,j{E) = M,dij - tij{E) , (2.10) 

where 6ij denotes the Kronecker delta. For this, we use the relations (2.7) to rewrite the 
integrand in (2.8) as follows: 

J2 [MEME) - fj(E)t,j{E)] = ME)M, - J2 fj(EME) 
j j 

= fj{E) {MAj - Uj{E)) . (2.11) 

Hence, one can rewrite the currents in the following form (i e {L, R, 1, . . . , N}): 

fj{E)T,j{E)dE, (2.12) 

1 

= hX^i, fjiE)'^ijiE)iE - fli) dE . (2.13) 

From (2.7) one immediately obtains the following properties: 

J2 ^iAE) = , Vj and Yl r^i(^) = , Vi . (2.14) 

i j 

Remark 2.5. Note that rij{E) ^ rji{E) in general. 

Using the relations (2.14) one immediately sees that the expressions (2.12)-(2.13) for 
the currents satisfy the conservation of charge and energy, respectively, i.e. 



/, = and ( ^^ + -h) = . (2.15) 



Remark 2.6. The second relation in (2.15) emphasises the fact that, in general, energy is 
conserved but not heat, in the sense that Jj 7^ in general. However, as we shall see in 
the next section (see Remark 3.1), heat is conserved in the linear response approximation. 
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2.2 The Strategy 

In Section 4, we will fix the values (Tl, /il) and (Tr, /ir), and determine the parameters 
(Tj, /ij), for i = 1 . . . , N, for which the self-consistency condition is satisfied: 

li = and Ji = for i = 1, . . . , iV . (2.16) 

Actually, this condition will be solved only in the linear response approximation. Then we 
will work out the remaining electric and heat currents: Jl = — /r and Jl = — Jr (these 
equalities are satisfied in the linear response regime; see Remark 3.1). Observe that in this 
situation, charge and heat dissipation occurs mostly at the left and right boundaries of the 
system. Indeed, due to the self-consistency condition, the dissipation occurs at first order 
at the boundaries while at most as second order effects along the system. 

Up to that point, the system is described by any scattering matrix S and has therefore 
no specific geometry (see Figure 4). In order to obtain a system with the linear geometry 
represented in Figure 3, and to interpret Jl and Jl as the currents flowing through the 
system, from left to right, we will consider that the global scattering matrix S is given in 
terms of the local (linearly ordered) scattering matrices, S'^^\ . . . , S'^^\ associated to the 
N quantum dots. 

In Section 5, we will study numerically the properties of the composite scattering matrix 
S, using random matrix theory (RMT), and show that the currents through the chain of 
quantum dots with self-consistent reservoirs are typically (in the sense of RMT) governed 
by Ohm and Fourier laws. 

3 General Properties 

In this section, we consider that the system is coupled to + 2 terminals (we use L and R 
to be consistent with the next sections) and suppose that its transport properties are given 
in terms of some scattering matrix S (see Figure 4). 




2 

Figure 4: The general system coupled to A^ + 2 terminals. 

In the linear response regime and under the assumption that the scattering matrix S 
does not depend on the energy, we will determine the general transport properties of the 
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system. All these properties are expected to hold on physical grounds and actually some 
of them have recently been obtained in full generality [37]. Nevertheless, in this paper, we 
are only interested in some particular situations in which the derivations of these properties 
are somehow simpler. 



3.1 Linear Transport Analysis 

In the remainder of this paper we shall always rely on the following. 
Assumption A3. We assume that, for i E {L, R, 1, . . . , N}, 

Ti = T + 5Ti, (3.1) 
/Xj = + , (3.2) 

where T > and // G M are some reference values, while 6Ti G M and dfii G IR are some 
"small" perturbations. For definiteness, we assume that all subsequent expressions are 
first order expansions in 5% and Sfj,i. 

To study the transport, we consider the first order expansion of any distribution func- 
tion: 

g{E; T„ /i,) = g{E; T, /x) + ^{E; T, /x) 6T, + -^{E; T, /i) % . (3.3) 
Note that the three considered distribution functions (2.1)-(2.3) are of the form: 

g(E;T^,,,.)^g(^E^^ . (3.4) 

Hence, one can rewrite the derivatives of g, with respect to Tj and Hj, in terms of the 
derivative of g with respect to E, and obtain the following result: 



df 

f^{E) = f{E; Tj, = f{E; T, ^x) - ^{E; T, ^x) 



E- fi 
T 



STj + 5iij 



(3.5) 



Substituting this first order expansion into the expressions (2.12)-(2.13) for the currents, 
and using the relations (2. 14), one sees that all terms containing f{E; T, /x) vanish and one 
is left with 

I, = ^L^^ + lJff^. (3.6) 



e T 
5T 



J, = ^4.)^^ + i<J)^, (3.7) 



where 



4°' = -jr§^(^-^'f'>^'i('^)'"^- P.8) 
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Remark 3.1. Using the relations (2.14), one sees that the first order approximations (3.7) 
for the heat currents satisfy 

J2J^ = 0. (3.11) 

i 

In Section 4, we will impose the self-consistency condition, li — Ji — for i — 1, N, 
and the conservation laws (2.15) and (3.11) will permit to interpret /l = — /r and Jl = 
— Jr as the electric and heat currents flowing through the system (from left to right). 



3.2 The Transport Matrix 



Here we establish some basic properties of the transport matrix. From the expressions 
(3.6)-(3.7) for the currents we are naturally led to the following definition. 



Definition 3.2. We denote by L the transport matrix: 



L(o) 



where L^'^^ = 



4? 



\-^RL 



T (fc) 

T(fe) 
^11 



(fe) 

-(fc) 

^Rl 



(fe) 
' LN 



(fe) 
■'NN 



r(fc) 
-^IR 



rik) 

^NR 
Ak) , 
^RR/ 



(k) 

A quick glance at the expressions (3.8)-(3.10) for the transport coefficients , re 



membering that Tij{E) = Midij — tij{E), shows that the Onsager relations hold, i.e. if 
S{E) is symmetric (for each E > 0), then L is symmetric. More generally, if the system 
satisfies the micro-reversibility property Sij^mniE, B) — Sji^nm{E, —B), where B is some 
applied magnetic field, then Lij(B) = Lji{—B). 

In general, the scattering matrix S depends on the energy. Nevertheless, in this pa- 
per, we shall restrict our attention to the energy-independent situations. The assumption 
that S does not depend on the energy is a good approximation in some limiting cases 
(e.g. looking at (3.8)-(3.10) one sees that for electrons at low temperature one may con- 
sider S{E) = S{Ep), where E-p is the Fermi energy, since {~df /dE){E) ^ S{E — Ep)) 
and leads to interesting consequences. In particular, we will see that the relations among 
the transport coefficients, L^^^ and L^f , will be very simple and will lead to some 
universal transport properties. For definiteness, we assume for the rest of this paper that 
the following holds. 



Assumption A4. The scattering matrix S does not depend on the energy. 

In order to simplify some derivations, without restricting our results too much, it is 
convenient to assume the following. 

Assumption A5. We assume that tij ^ Ofor all i, j E {L, R, 1, . . . , N}. 

Remark 3.3. Recalling the properties (2.7), one sees that the Assumption A5 implies that 
tij 7^ min{Mi, Mj}, for alH, j G {L, R, 1, . . . , A^}, and this permits to avoid total back- 
scattering (i.e. tii 7^ Mi for all i). Although not optimal, the Assumption A5 also insures 
the existence and unicity of the solution of the self -consistently condition (2.16). 
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Under the assumption that the scattering matrix S does not depend on the energy, one 
sees in the expressions for the transport coefficients (3.8)-(3.10) that some integrals of the 
following form appear (with n = 0, 1, 2): 

More explicitly, 

C^^{n)= x'^e--" dx and C^(n) = -:r dx , (3.13) 

where the signs + and — correspond to the Fermi-Dirac and Bose-Einstein cases, respec- 
tively, and xq = —ji/ {k^T) is some reference parameter. Here, e K in MB and FD, but 
Xo > in BE. 

Remark 3.4. Although C{n) depends on xq, we will see that the transport properties of the 
system are essentially independent of xq. This explains why we do not write C(n, xq). 

In terms of C{n), one can write 

Llf = ^C(0)r,,-, (3.14) 
Lg) - ^fcBTC(l)r,,, (3.15) 

Lg) = i^C(2)r.,. (3.16) 

In Appendix A, we show that C(0), C(l) and C(2) are positive for all 

^-l,oo)inMB 

Xo = --4;-e<! (-CX), oo) in FD . (3.17) 



(0, oo) in BE 



In the Maxwell-Boltzmann case, when xq e (-oo, -1], we find that C^^(O) and C^^^(2) 
are positive, but C^^^(l) is non-positive. As one may easily check, all the subsequent 
results also hold in this situation. Nevertheless, in order to be able to treat the three consid- 
ered cases on the same footing we make the following restricting assumption. 

Assumption A6. In the Maxwell-Boltzmann case: Xo = —ii/{k-BT) G (—1, oo). 

Remark 3.5. Note that the coefficient C^^(O) diverges as xq — > 0. This may be interpreted 
as a Bose-Einstein condensation, and shows that our theory breaks down in describing 
properly the transport of Bose-Einstein condensates. 

Remark 3.6. In mesoscopic physics, one is usually interested in electronic transport at 
low temperature, i.e. xq = —iJ,/{kBT) — > — oo. Observe that in this limit C^°(l) 
and consequently lJ^^^^ — > 0. In order to have non-zero transport coefficients L^f^^^ in 
this regime, one usually considers the first order term in the Taylor expansion, Tij{E) — 
+ ^.{^,){E - /.) + 0{{E - which leads to L™^^) - ^ .(/.). 
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Remark 3.7. Observe that the transport coefficients L-^^ depend on the properties of the 
quantum dots through the scattering matrix S. Here we can see why the assumption A5 
is useful. Indeed, suppose that tu = Mi for some i (total reflection). Then tij = for 
all j i and consequently — for all j. In other words, in such a situation the i-th 
reservoir is disconnected from the rest of the system. 

Looking at the relations (3. 14)-(3. 16), one sees that the transport coefficients are related 
as follows: 



I^-T^m,^ I^-Q^I^ and L^^Q.L^, (3.18) 



h 

where 

«. = ^§|i|>0 ana «. = ^§||>0. C3..) 

Remark 3.8. The quantity Qi is related to the Seebeck coefficient (see Remark 4.3). In the 
Fermi-Dirac case, Mmx^^^oo Q2/T'^ = (7r^/3)(A;B/e)^ is called the Lorentz number. 

Recalling that Tij — Mi5ij — Uj, one immediately obtains the following properties 

(A; = 0,1, 2): 

and 

J2 Lf = , Vi and ^ if = , Vi . (3.21) 

i 3 

As we shall see later on, the following ratio, IZ, will play an important role. Using 
standard techniques, we show in Appendix B that (for all xq): 

n=-%^ , ^^^^ e (0, 1) . (3.22) 
VQ^ VC(0) • C(2) ^ ' ^ 



\C{1) C{2)J' ^^-^^ 



Remark 3.9. In terms of the matrix 

C 

the inequality < 7^ < 1 is equivalent to det(C) = C(0)C(2) - C{iy > 0. 
3.3 The Entropy Production 

In this subsection, we show that the transport matrix L is real positive semi-definite. In 
other words, we show that the entropy production rate is non-negative: 

2N + 4: 

as = ^^^^^^ ^ ' (3.24) 
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where the thermodynamic forces are 

( S/ii^i/e if i^l,...,N + 2 

\STi^^N+3)/T if i = TV + 3,..., 27V + 4 " ^^'^^^ 

Here, we have set L = and R = iV + 1. Let Xi = 5ni/e and Zi = STi/T, for 
i e {0,1, . . . , N, N + 1}. Then, in Appendix C, we show that one can rewrite the entropy 
production rate as follows: 

N+l 



E (3.26) 



i,j = 
i <j 

where 

lij = — XjY + {Zi — ZjY — 2TZCij , Cij = XiZj + XjZi — XiZi — XjZj . 

Observe now that {—hf^) > for all i j. Therefore, to show that > it is sufficient 
to show that 1^ > for all i 7^ j. As we shall see, it will be crucial that the ratio TZ satisfies 
< 7^ < 1. Assume first that Qj < 0, then (TZ > 0) 

> {Xi - Xj)^ + {Zi - Zj)^ > . (3.27) 

Assume next that Cij > 0, then (TZ < 1) 

lij > {Xi - Xj^ + {Zi - Zj^ - 2Cij = {Xi - Xj + Zi- Zjf > . (3.28) 

This ends the proof that as > 0. 

Remark 3.10. The term (Xj — Xj^ accounts for the entropy production due to the electric 
current, and it is well known in mesoscopic physics (see e.g. [35]). The term {Zi — Zj^ is 
related to the heat current, and one sees that there is also a thermoelectric term, —2TZCij, 
which might take positive and negative values. 



3.4 Equilibrium and Non-Equilibrium States 

Here we present some equivalent characterizations of the equilibrium and non-equilibrium 
states of the multi-terminal system. We start with the following. 

Definition 3.11. We say that the system is at equilibrium if 

7l = Ti = • • • = Tat = Tr and //l = //i = • • • = //^ = //r . 
Otherwise, we say that the system is out of equilibrium. 

Now, as one expects, the system is at equilibrium only in the situations in which all the 
electric and heat currents vanish, i.e. 

{System is at equilibrium} <^=^ {/^ = and Jj = 0, Vi} . 
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This equivalence is obtained by writing the no-current condition, Jj = Jj = for all i e 
{L, R, 1, . . . , A^}, in a matrix form and by using the properties of the transport coefficients 

(k) 

Lij as well as the fact that <7l < 1. The details are presented in Appendix D. 

Also, one can characterize the equilibrium and non-equilibrium states of the system in 
terms of the entropy production: 

{System is at 1 a — and / System is 

equilibrium j ^ \ out of equilibrium 

These equivalences easily follow from the expression (3.26) for ag and the inequality 
(3.27). The details can be found in Appendix D. 




(J, > . 



4 The Self-Consistency Condition 

We now turn to the resolution of the self-consistency condition (2.16). Using the expres- 
sions (3.6)-(3.7) one can rewrite the self-consistency condition as follows (i = 1, . . . , iV): 

t (4»' + !<]• f ) = - (l<»> ^ + !<]• f ) , (4.29) 
t (Uff '-f + US' ^) = - E (4' ^ + if ^) ■ (4.30) 

j=l ^ ^ j=L,R ^ ^ 

We recall that we are given the values (Tl, /il) and (Tr, yUa), so that the right-hand side 
of the above equations are supposed to be known (Tj = T + 6Tj and fxj = n + Sjij). 
Therefore, the self-consistency condition constitutes a set of 2N equations for the 2N 
unknown variables Ti, . . . , Tat and fxi, . . . , hn. 

To solve these equations, it is convenient to introduce the vectors X,Y & and the 
N X N matrix defined by = 1,..., N): 

X, = ^, y, = ^, and iL^h = L^. (4.31) 
Remark 4.1. Recalling that is the (iV + 2) x (A^ + 2) matrix defined in Definition 3.2, 

^(0 

We also define for £ — L,R 



one sees that is the reduced matrix associated to the "central" terminals 1, . . . , A^. 



Xe^ — , = — , and De^ 
e 1 



(4.32) 



Then the equations (4.29)-(4.30) can be rewritten as follows: 

lP {X + Q,Y) = -J2{Xe + QiYe) De , (4.33) 



^=L,R 



Lg) {Q,X + Q2Y) ^ Yl (^1^^ + ^2>^^) De . (4.34) 



e=L,R 
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Since TZ ^ 1, one has Q2 7^ (Qi)^- As a consequence, if one takes (Q2/Q1) times (4.33) 
and subtract (4.34) and similarly if one takes (— l/Qi) times (4.34) and add (4.33), one 
obtains 

L^X DeXe and L^V = - J] DeYe . (4.35) 

We see here that the equations (4.33)-(4.34) decouple and that the equations for the chem- 
ical potentials and for the temperatures are actually identical. In order to proceed, we need 
to know that the matrix is invertible. This fact is proved in Appendix E by showing 
that is real positive definite. In particular, this shows that the self-consistency condition 
(2.16) has a unique solution. 

By inverting the matrix Lq\ we obtain the following self-consistent temperatures and 
chemical potentials (i — 1, . . . , N): 

T, = rL + A(7R-rL), (4.36) 

/^j = /^L + (//R - //l) , (4.37) 

where 

N 

A = 5](rc')i,- , with rc = (r,,)5=i. (4.38) 

i=i 

The details about their derivation are given in Appendix F. 

By substituting these expressions into the relations (3.6)-(3.7) for the currents, with 
i — L, and by using the properties of the transport coefficients Llj', one deduces the 
(average) currents flowing through the system (from left to right): 

Il = 0-0 ( ) + f — ) ' (4-3^) 

Jl = a,[^^^^]+aJ^^] , (4.40) 



where 

N 

(7o = 4r + J]Aj4°- , (Ti = Qi(To and (72 = Q2f7o. (4.41) 

Remark 4.2. The corresponding expressions to (4.36)-(4.41) in the case of electronic trans- 
port at low temperature (Tj = and Ji — for i e {L, R, 1, . . . N}) were obtained in [27]. 

Remark 4.3. Setting Jl = 0, one deduces the thermoelectric field £ = —W = iSVT, 
where S = Qi/T = ks/e ■ C{1)/C{0) > Ois the thermopower or Seebeck coefficient. 
Note that S does not depend on N. For further discussions see e.g. [39-45]. 

Using the relations (4.36)-(4.37), one may rewrite the currents (4.39)-(4.40) as follows: 



16 



Philippe A. Jacquet 



Although the multi-terminal system considered so far has no specific geometry (see 
Figure 4), we shall nevertheless use the word "profile" for the arrangements ni, . . . , hn 
and Ti, . . . , T^. This terminology will be fully justified in Section 5, where we will restrict 
our attention to the linear chain of quantum dots (Figure 3). 

Remark 4.4. Note that the temperature and chemical potential profiles are decoupled, i.e. Tj 
only depends on (Tl,Tr) and jii only depends on (/iL,AiR)- On the other hand, since 
Qi > 0, there can be a heat current due to a gradient of chemical potential (Tl = Tr), and 
reciprocally, there can be an electric current due to a gradient of temperature (//l = A*r)- 

Remark 4.5. We see that the form of the profiles is given by the coefficients Ai, . . . , An- 
Hence, in what follows we will use the word "profile" to refer to these coefficients. Looking 

at the expression (4.38) for Ai one sees that it does not depend on the distribution function 
(2.1)-(2.3) describing the reservoirs. In this sense, we say that the profile is universal. One 
may wonder whether this is true for any distribution function /. It seems not to be the case. 
Indeed, the key ingredients to obtain the preceding results are: (i) the specific form (3.4) 
of / and (ii) the ratio TZ, defined in (3.22), must satisfy < 7?. < 1. Basically, we expect 
that any distribution function / satisfying (i)-(ii) will lead to the same results. Note that 
universal temperature profiles were also found for a quantum harmonic chain coupled at 
both ends to two phonon reservoirs at different temperatures [46]. 

Remark 4.6. Although the profiles are the same in the three considered situations, one can 
distinguish between the classical and quantum chain of quantum dots. Indeed, assume we 
are given the scattering matrices, S'^^\ . . . , S^^\ associated to the N quantum dots. Then, 
as explained in detail in Appendix H, one can compose them into a global scattering matrix 
S, out of which one extracts the transmission probabilities tij (see (2.6)). This is the natural 
quantum way of working. 

In order to extract the interference effects, one may first compute the probability ma- 
trices, P''^\ . . . ,P'^^\ defined by Pi^l^j^ = |5'^j;LriP' then compose them into a global 

probability matrix P and finally set tij = ^^li ^n=i Pij;mn- This way of computing tij 
will be referred to as classical. 

As one may check, the composition law of the probability matrices is the same as for 
the scattering matrices, i.e. all expressions in Appendix H hold with P instead of S. In this 
sense, we will study in Section 5 the classical versus quantum situations (see also [47,48]). 

In Appendix G, we derive an interesting alternative expression for the coefficients A^. 
We obtain (i = 1, . . . , A^) 



where rc(j, i) denotes the (j, i) minor of Fq. 

Remark 4.7. As an illustration, let us apply the relation (4.44) in the case N — 1. One has 



E;=i(-l)^+^det(rc(j,i)) 



(4.44) 



Ef=i(-l)*+^ det(rc(j,0) [t,L + t,R] ' 




t 



tvL IJ-h + ^IR /J-R 



(4.45) 
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This relation for fxi corresponds to the one obtained in [35] for a three-terminal conductor 
under the conditions Tl = Ti = Tr = (which implies Ji = 0, since 6Tj = and 
L™^^^ = in this case) and Ji = 0. 

An important consequence of (4.44) is (i — 1, ... , N): 



This implies that the temperatures and chemical potentials of the reservoirs, connected self- 
consistently to the system, are bounded between the boundary values Tl, Tr and /^l, A*r> 
respectively. In other words, for i = 1, . . . ,N, one has 

min{TL, Tr} < Tj < max{TL, Tr} and min{/XL, A*r} < A«i < max{//L, A^r} ■ 

Remark 4.S. The property (4.46) also implies that the self-consistent parameters Ti, . . . , T^v 
and //i, . . . , //AT are increasing functions of Tl, Tr, and //l, I^r, respectively. 

Moreover, looking at the expressions (4.41) for the coefficients ctq, cti and a2, and 
recalling that L^f < if i 7^ j, Qi > and Q2 > 0, one deduces that for A; = 0, 1, 2: 



Remark 4.9. There are some simple consequences of (4.47) with respect to (4.39)-(4.40). 
(i) If /^L = A*R, then the heat current flows from the hot reservoir to the cold one, and 

similarly (ii) if Tl = Tr, then the electric current goes from the high chemical potential to 
the lower one. Note, however, that by setting ^Lj /Wr and Tl, Tr to appropriate values there 
might be a heat current going from the cold reservoir to the hot one. (A similar statement 
holds for the electric current.) 

4.1 Ohm and Fourier Laws 

From the relations (4.39)-(4.40) for the electric and heat currents, one can find explicit 
expressions for the electric and heat conductances. If Tl = Tr, then setting /i = eV, 
where V denotes the electric potential, one obtains 



A, e (0, 1) . 



(4.46) 



(jfe < . 



(4.47) 



Ge (^R - Vl) , 



(4.48) 



where the electric conductance Ge is given by 



Ge — —(To . 



(4.49) 



If /l = 0, one deduces 



Jh — —Gh (Tr — Tl) , 



(4.50) 



where the heat conductance Gh is given by 




{Q2 - Ql) 



(Jo = 



fc|T G(0)C'(2) -G(l)^ 

G(0)2 



(7o . 



(4.51) 



croT 



T 
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Recalling that < 7^ < 1, one sees that Q2 > {QiY- Consequently, since (Tq < 0, one 
deduces that the electric and heat conductances are positive: 

G'e>0 and Gh > . (4.52) 

Assume now that the system has the linear geometry represented in Figure 3 and let a 
denote the length of the quantum dots and S denote the width of the openings connecting 
any two dots (i.e. the effective cross section of the system). Since L = Na corresponds to 
the length of the system made of N quantum dots, one may introduce the global gradients: 

VT/ = _R___t and VT = . (4.53) 

ij ij 

Then, under the same conditions as in (4.48) and (4.50), the current densities, Xl = /l/S 
and = </l/S, are given by the global Ohm and Fourier laws: 

2l = -K^W and Jl = -^hVT , (4.54) 

where the electric and heat conductivities, and k\^, are given by 

LGp LGh 
= ^ > and = — ^ > . (4.55) 

Remark 4.10. In the Fermi-Dirac situation, one has 

limQi^Oand limQ^^-f— =^ lim '^'^ _ ^V^b 



a;o-»0 lo^O 3 \ e / a;o->0 n^T 3 \ 6 

The last relation is the Wiedemann-Franz law giving the Lorentz number. 
In a similar manner, one may introduce the following local gradients: 

VV{i) = ~ and VT{i) = ~ . (4.56) 
a a 

Let X{i) = Xl and J'(i) = Ji^, where i E {1, . . . , iV — 1}, denote the net local current 

densities flowing from the i-th quantum dot into the {i + l)-th quantum dot. Then, using 

the equations (4.36)-(4.37) and (4.54), one obtains the local Ohm and Fourier laws: 

= _fi;^(i)v\/(i) and J{i) = -fi;h(i)VT(i) , (4.57) 

where the local electric and heat conductivities are given by 

Remark 4.11. Although the global conductivities, and k^, are always positive, we will 
see, in Section 5, that the local conductivities, Ke{i) and may take negative values. 

The coefficients ao, ai and a2, the conductances Gc and Gh, and the conductivities 
and K,h depend on the distribution function / describing the reservoirs. Nevertheless, they 
are all proportional to ctq, with some multiplicative factor independent of A^. The important 
observation is that the dependence on the nature of the particles is entirely contained in 
these multiplicative factors. This permits to introduce the following universal quantity. 
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Definition 4.12. We define the universal conductivity as 



h 



N N 



k{N) 



e2C(0) 



(4.59) 



i=i j=i 



Remark 4.13. As an illustration, let us consider the case = 1. One has 



«(1) =iLR + 



(4.60) 



^IL + ^IR 



As expected, the universal conductivity is the sum of two contributions: a direct coherent 
transmission and an indirect incoherent transmission through the intermediate terminal. 
This result is well known in mesoscopic physics [35], and one sees that the relation (4.59) 
is the generalisation to N self-consistent reservoirs. 

Looking at the expression (4.59) for the universal conductivity, one sees that k{N) > 0, 
since ctq < and C(0) > 0. Using moreover the relations (2.7) and ^lr < min{ML, Mr}, 
one finds the following optimal bounds: 



5 Numerical Simulations 
5.1 The Chain of Quantum Dots 

To prove the validity of the global Ohm and Fourier laws in our model, it remains to show 
that the transport is normal, that is, when the system becomes macroscopic the global 
electric and heat conductivities, He and k,-^, become intrinsic properties of the macroscopic 
system and consequently become independent of its length L — Na. More precisely, 
one has to prove that the global electric and heat conductivities converge to finite limits 
as L — > oo. Looking at the relations (4.55) and (4.59), one sees that and are both 
proportional to k{N), with some multiplicative factor independent of N, and consequently 
it is sufficient to investigate the limit of k,{N) as A?^ — > oo. 

As explained previously, the multi-terminal system considered so far has no specific 
geometry (Figure 4), and therefore k{N) will not converge as ^ cxo in general. In 
order obtain a finite limit for k{N) as ^ oo we shall proceed as follows. 

Let (5'*^''^)^^^ be an infinite sequence of scattering matrices, each one being associated 
to a quantum dot, and let Sn denote the composite global scattering matrix associated to 
{S^^\ . . . , S^^^) . The details concerning the construction of (S'jv are given in Appendix H. 
Having set-up the linear geometry represented in Figure 3, one may expect that k{N) will 
converge to a finite value, as A^ ^ cxo, and this is what we found in typical numerical 
experiments (see below). Unfortunately, the composite scattering matrix, S^, is given in 
a rather cumbersome form in terms of 5"^^^ , . . . , S^^\ which prevents us from proving the 
existence of hmjv-*oo «(A^). 



< k{N) < N ■ min{ML, Mr} . 



(4.61) 
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Nevertheless, we will see by using numerical simulations under random matrix theory, 
that the universal conductivity k{N) typically admits a finite limit, {k°°), as iV — > cxd. This 
shows that the global Ohm and Fourier laws hold on statistical average, and we may write 

(c) = I ^ in (5.1) 

and 

_ a klT C{0)C{2)-C{ir 

Moreover, we will see that typically the profile, Ai, . . . ,An, becomes linear as the 
number of dots N ^ oo. Introducing the variable x — i/{N + 1) e [0,1], this means that 
in the thermodynamic limit N ^ oo, one has 

{A{x))^x, (5.3) 

and consequently 

{Vix))^Vj^ + iVn-Vj^)x and {T{x)) + (T^-Tj^) x . (5.4) 

Hence, the local Ohm and Fourier laws also hold on statistical average, with position- 
independent electric and heat conductivities. 

Remark 5.1. Here we compare the values of the average electric and heat conductivities for 
the three considered situations. Indeed, looking at the expressions (5. l)-(5.2), one sees that 
they only differ from a multiplicative factor, which depends on the nature of the particles. 
In the electric and heat cases, we found the following order: BE > MB > FD, for all xq > 
(the common admissible region). 



5.2 Random Matrix Theory 

As mentioned in the introduction, the classical EY-model describing particle and energy 
transport (Figure 1) has the following important property: it is chaotic. More precisely, 
the cells, without the discs, which are making up the system are ergodic. In analogy, we 
assume the following. 

Assumption A7. The quantum dots are classically chaotic. 

Under this assumption, it is natural to describe the transport properties of the quantum 
dots with random matrix theory (RMT) [49]. More precisely, we will consider that the 
quantum dots are described by the Wigner-Dyson circular ensembles (with /3 = 1,2). 

Here, we restrict our attention to the one channel situation, i.e. we assume that Mi — 1 
for z e {L, R, 1, . . . , N}. Now we consider that the 3 x 3 complex random matrices, 
S^^\ . . . , S^^\ are independent and identically distributed over U (3) with some measure. 
More specifically, we consider the following well-known cases: 

(1) The Circular Orthogonal Ensemble (COE, (3 — 1) describes the case of time-reversible 
systems. 
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(2) The Circular Unitary Ensemble (CUE, jS = 2) corresponds to the case in which the 
time-reversal symmetry is broken. Such a situation appears, for example, when the 
system is submitted to some external magnetic field. 

In the sequel, we shall be interested mainly in the following quantities: 

(i) The universal profile Ai, . . . , Apf, which gives the shape of the temperature and 
chemical potential profiles. In the subsequent plots, we naturally set = = 
and An+1 = = 1. 

(ii) The universal conductivity k,{N), which gives the electric and heat currents through 
the chain of N quantum dots. 

As explained previously, we shall compare the classical and quantum situations by working 
in terms of the probability matrices, P^^^ , . . . , P'^^\ and scattering matrices, S'^^^ , . . . , S'^^^ 
respectively (see Remark 4.6). 

For a review on RMT we refer to [49]. To implement the numerical simulations, we 
have followed [50]. In the sequel, the averages are made over an ensemble of lO'^-lO^ 
realizations, depending on the size N of the system. 

5.3 The Global Transmission Probabilities 

We denote by (^f )i and (^f )2 the average of the transmission probabilities in COE 
(P — 1) and CUE (/3 = 2), respectively. It turns out that one can determine their exact 
values [49]: 



Now we use the construction explained in Appendix H to build the scattering matrix S 
associated to quantum dots. As the reader may guess, it is then very difficult to make any 
analytical statement concerning the averages (tij)i and {tij)2- We thus turned to numerical 
simulations, and made the following observations for large system sizes N: 

(1) The classical and quantum global transmission probabilities, tf^ and tf^, are dijf ev- 
ent, but their averages, {tf^p and {t^^)i3, are the same, showing that on average the 
interferences are negligible when N is large. 

(2) The averages {tij)i3 do not depend on the system size N. 

(3) They are symmetric: (tij)/3 = (tji)/3. 

(4) The average couplings between the terminals, (t^)^, depend on\i — j\ and are short 
ranged, i.e. {tij) ^ ~ if |i — j | > 2. In particular, the average probability matrix has 




(5.5) 
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essentially the following form: 



* * * 











* * 

* * 



The values of the matrix elements represented by are smaller than 0.05. 

A typical example illustrating the above observations is given in Figure 5, where = 
20. In this figure, five terminals are considered: i E {L, 5, 10, 15, R} (with L = and 
R = 21). The different curves correspond to {tie)/3 for some i E {£ — 4, ...,£ + 4}. For 
example, the middle curves give the following values: (t7,io)/3, • • • , {tio,io)/3, • • • > (^i3,io)/3- 




Figure 5: The (classical and quantum) average transmission probabilities, {tii)/3, for i E 
{L, 5, 10, 15, R} (with L = and R = 21) and some i G {£ - 4, ...,£ + 4}. For example, 
the middle curves give: (t7,io)/3, • • • , (^io,io)/3, • • • , (^i3,io)/3- Left: COE (/3 = 1). Right: 
CUE (P = 2). 



5.4 The Profiles and Currents 

Let us recall the expressions for the coefficient and for the universal conductivity k,{N) : 



N 



A = - J2^rc% and k{N) = N 



N N 



In Figure 6, we show the average profiles {A{x))p for A = 5. If Tl < Tr, these curves 
correspond essentially to the temperature profiles along the system. (The same holds for 
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the chemical potential profiles if /xl < AfR)- We see that all these profiles are not linear, 
the non-linearity being more stressed in the quantum situation, showing the effect of the 
interferences. However, we found that all these profiles, while conserving their shape, get 
closer and closer to linear as the number of dots increases and that eventually they become 
linear in the limit N ^ oo. 




Cell X = i/(N+l) Cell x = i/(N+l) 

Figure 6: The classical and quantum average profiles {A{x))f3 for = 5. The diagonal 
lines are added to guide the eye. Left: COE ((3 = 1). Right: CUE (/3 = 2). 



Remark 5.2. In order to emphasise the role of interferences, let us consider the following 
simple toy models. We consider that the classical probability matrices are identical and 
given by Pj;^ = 1/3, for all i, j. In such a situation, the profile A{x) is linear (for all N). 
On the other hand, the following unitary scattering matrices 

1 6 3 e 3 \ 
e^^ 1 M , (5.6) 
e^^ ^ 1 / 

which satisfy the equiprobability property, |5'fj^^p = |, do not lead to a linear profile, 
showing the interference effects (see Figure 7). A discussion about the conductivities in 
these cases is given in Section 5.5. 



Remark 5.3. One may have noticed that the average temperature profiles differ from those 
in the EY-model. Indeed, in the EY-model the temperature profiles were either linear, 
convex or concave functions of the cell position, depending on the values of the injection 
rates 7l and jr. Basically, the reason was that, in the presence of a particle current, for 
example, going from right to left, the disc i hears more often from cell i + 1 than from 
cell i — 1. It has therefore a greater tendency to equilibrate with the right than with the 
left, causing a curvature in the temperature profile. In particular, the temperature profile 
depends on all the boundary values Tl, Tr, 7l and 7r. In the present model, the coupling 
among the terminals is governed by the scattering matrix S and we have assumed it is 
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Cell X = i/(N+l) 

Figure 7: The profiles of the toy models in Remark 5.2 (N — 5). 

energy-independent leading to a self-consistent temperature profile dependent only on S 
and on the boundary values Tl and Tr. However, we expect that in the general energy- 
dependent S-matrix situation, the temperature profile will also depend on /xl and hr. To 
see whether in this case the EY-model and our model agree needs further investigations. 

We have shown in the preceding section that, if /^l = A*r and Tl < Tr, then the heat 
current Jl goes from the hot reservoir (Right) to the cold one (Left), as one expects from 
the second law of thermodynamics. In Figure 6, one sees that the profiles are monotically 
increasing from left to right. This means that, on average, heat also flows locally from hot 
to cold, i.e. on statistical average, the second law of thermodynamics also holds locally. 

All the preceding results concern statistical average values. In Figure 8, we present 
different realizations and one sees that very different profiles may occur. Roughly speak- 
ing, by choosing an appropriate combination of the S-matrices, S^^\ . . . , S^^\ one may 
generate almost any profile. 

In particular, one sees in Figure 8 that there exist particular realizations for which the 
heat current flows locally from cold to hot. In other words, the profiles are not monotone 
increasing for every realization, but only on statistical average. Such phenomena are well 
known in mesoscopic systems (see e.g. [35,51,52]). 

We observed that the quantum dots about which such phenomena occur (i.e. such that 
Ai < Ai_i and consequently Ke(i — 1) < and K\^{i — 1) <{)) have typically a very high 
reflection probability tu. As the system size N grows, we see that such phenomena dimin- 
ish in intensity and frequency. Basically, we expect that, when the system is macroscopic, 
the second law of thermodynamics also holds locally in every realization. 

Moreover, in the thermodynamic limit N — oo, we see that the average profiles be- 
come linear in all situations (COE, CUE, classical and quantum): 

{A{x))p^x, Vxe[0,l]. (5.7) 

Consequently, 

{V{x))p^V^ + {V^-V^)x and {T{x))p^T^+ {T^-T^) x . (5.8) 
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Cell x = i/(N+l) 
Figure 8: The quantum profiles for some realizations (A^ = 5). 



We next present some results concerning the universal conductivity k,{N). In Figure 9, 
we show {K{N))fs for small N. Observe that in the COE cases the conductivity decreases 
with N while in the CUE cases it first increases and only then decreases. 




^ CUE Classical 
^ CUE Quantum 

COE Classical 
COE Quantum 



Figure 9: The average universal conductivities {k{N))i3 for small N. 
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For larger values of A^, we see in Figure 10 that the average universal conductivity 
seems to admit a finite limit {K°°)f3 as iV — > oo. To support this observation and obtain 
a numerical value for {K,°°)f3, we have proceeded as follows. We found that the curves in 
Figure 10 are very well fitted (the largest chi-square being 0.035) by functions (not shown 
in the figures) of the following form: + Cp/N"'. Here, cg is some positive constant 

and a is an exponent found to be the same in COE and CUE, and given by 0;*^' — 1/2 and 
CK^" = 1. Moreover, 
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Figure 10: The average universal conductivities {k,{N))p as a function of N. 



Let us denote by and {K^)f3 the corresponding limits: 



{'^h //? - ^ T IP ■ (p.W) 



and 

_a 4T C(0)C(2) - C(1)^ 
E ~li C(0) 

These numerical results show that the global Ohm and Fourier laws hold on statistical av- 
erage in our chain of quantum dots with self-consistent reservoirs. Since the average tem- 
perature and chemical potential profiles become linear in the limit ^ cxo, it follows that 
the local Ohm and Fourier laws also hold on statistical average, with position-independent 
electric and heat conductivities, i.e. 

{I{x))p^-{K^{x))pV{V{x))p and {J{x)) ^ ^ -{k^{x)) {T{x)) p , 

with {Kf{x))ii = {K^)f3 and {n^{x))fs = {k^)i3 for all x e (0, 1). Here we have kept the 
global gradients (4.53) constant while taking the limit N ^ 00. 
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Remark 5.4. In the single channel case (Mj = 1, for all i), the optimal bounds (4.61) 
become < k{N) < N. Let us consider the two extreme situations in which all the 
quantum dots are identical. 

(1) Total back-reflection: 



This corresponds to the case in which all the intermediate terminals are disconnected 
from the system and the transmission through the system is optimal (i.e. ^lr = ^rl = 
tn — • • ■ — tj^N — 1 and all other — 0). We found k{N) ~ N, for all N. Hence, 
a perfect conductor has infinite electric and heat conductivities in the limit — > oo. 

Although we have shown that typically (in the sense of RMT) k{N) converges to a 
finite limit as A'^ — > oo, one deduces from the examples in Remark 5.4 that the linear 
geometry, represented in Figure 3, is not sufficient to ensure the existence of such a limit. 
The characterization of the sequences of scattering matrices, (5'*^'^^) for which the limit, 
limAT^oo f^iN), exists and is finite remains an open problem. Note however that, approxi- 
mate results were obtained in some particular cases [23,27]. Note also that the particular 
ordered cases, in which all the local scattering matrices are identical, are not simpler to 
handle with the relations given in Appendix H. 

In Figure 10, one sees two manifestations of weak localization: 

(1) The classical conductivity is larger than the quantum conductivity. Hence, the effect 
of the interferences is to decrease the intensity of the electric and heat currents. Note 
that this effect is rather subtle. Indeed, one has 



Since the average global transmission probabilities are the same in the classical and 
quantum situations, the effect of the interferences on the conductivity must be due to 
correlations among the random variables tij. 

(2) The application of an external field (CUE) increases the values of the conductivities 
and consequently increases the intensity of the electric and heat currents. This can be 
understood as follows: In equation (5.5) and in Figure 5, one sees that the presence 
of a magnetic field will (on average) increase the transmissions probabilities. 




In this case, we found k{N) 



for all N. 



(2) Optimal transmission: 





(5.11) 
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Remark 5.5. The case of electronic transport at low temperature was already considered by 
D'Amato and Pastawski [27]. Modelling the system as a nearest-neighbour tight-binding 
Hamiltonian they found approximate expressions for the global transmission probabilities 
tij, which turn out to satisfy the properties (2)-(4) presented in Section 5.3. From these 
approximate t^j they deduced that, in the limit the system size ^ oo, the self-consistent 
chemical potential profile is linear and the universal conductivity is finite, in agreement 
with our results. However, in our opinion, these (nice) results need further analysis to 
constitute a rigorous derivation of Ohm's law in their model. 

5.5 Chaos, Disorder and Decoherence 

In this paper, we were not interested in finding the optimal conditions to obtain Ohm and 
Fourier laws, but rather to establish a quantum version of the classical chaotic EY-model. 
This led us naturally to use RMT to characterize the quantum dots. Nevertheless, one 
may wonder whether chaoticity or randomness, as modelled by RMT, is essential to obtain 
normal transport, i.e. for Ohm and Fourier laws to hold. The answer is no, since, for 
example, the two toy models presented in Remark 5.2, which may be interpreted as ordered 
non-chaotic cases, do lead to Ohm and Fourier laws. Indeed, one sees in Figure 1 1 that 
their corresponding universal conductivities converge to finite values as the number of 
dots increases. Observe however that, contrary to the RMT situations, the conductivities 
increase with the number of dots and the "quantum" universal conductivity is higher than 
the "classical" one. This shows in particular that, although interferences are typically (in 
the sense of RMT) destructive, they may be constructive in some particular cases. Actually, 
this phenomenon is typical in ordered chains (see below). 



4 




500 



N 



Figure 1 1 : The universal conductivities k{N) of the toy models in Remark 5.2. 
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Another natural investigation concerns the effects of disorder. We say that the chain 
is ordered if all the local scattering matrices are identical (i.e. S^^^ = ■ ■ ■ = S^^^) and 
disordered otherwise. The results presented in Figures 6 and 10 concern the disordered 
chaotic cases and in Figures 7 and 1 1 one can find two examples of ordered non-chaotic 
cases. The ordered chaotic situation is obtained by using RMT with the condition that for 
each realization S^^^ = ■ ■ ■ = S^^\ The results are given in Figures 12 and 13. Observe 
that by applying a magnetic field the "s- shapes" of the profiles are reversed and that the 
conductivities behave qualitatively as in the (ordered) toy model cases. Note also that in 
the quantum COE situation, the asymptotic value of the universal conductivity is about 
three order of magnitude larger in the ordered case than in the disorder one. This shows 
that disorder may have a very strong localization effect. 




Figure 12: The average profiles {A{x))i3 in the ordered cases (N = 5). The diagonal lines 
are added to guide the eye. Left: COE ((3 = 1). Right: CUE ((3 = 2). 




CUE Quantum 
COE Quantum 
CUE Classical 
COE Classical 



Figure 13: The average universal conductivities {k:{N))i3 in the ordered chains. 
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Finally, let us discuss the effects of decoherence. The idea is to view the self-consistent 
reservoirs as an effective environment acting on a chain of coherent quantum dots. For 
simplicity, let us assume that all the self -consistent reservoirs are equally coupled to the 
dots, i.e. 1 — |5'22V = A e [0, 1] for all k. (In the simulations, the values of the coupling 
parameter A are taken in small intervals of width 0.01.) Then, A = corresponds to the 
completely decoupled situation and A = 1 corresponds to the maximally coupled one. 

For any fixed value of A, we found that the profiles and conductivities are qualitatively 
similar to those obtained previously. Interestingly, when A increases, we observed that the 
conductivity increases in the disordered cases while it decreases in the ordered ones. As an 
illustration, we show in Figure 14 the results for three different values of A in the quantum 
disordered chaotic chain. In this case, the profiles (not shown) are similar to those obtained 
in Figure 6 and tend to get closer and closer to linear as the coupling parameter A increases. 




Figure 14: The average universal conductivity {k{N))i in the quantum disordered chain, 
for three different values of the coupling parameter A. 



6 Concluding Remarks 

We have presented a model for charge and heat transport using the Landauer-Biittiker scat- 
tering approach. In the linear response regime, we have seen that the Onsager relations 
hold. Then, assuming that the scattering matrix S does not depend on the energy, we have 
shown that the transport matrix L is real positive semi-definite and we have characterized 
the equilibrium and non-equilibrium states of the system in terms of the currents and in 
terms of the entropy production. In particular, we have shown that our multi-terminal sys- 
tem satisfies the first and second laws of thermodynamics. 
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In order to obtain an effective quantum version of the Eckmann- Young model with 
rotating discs, we have imposed the self-consistency condition: /j = Jj = for i = 
1, . . . , iV. This condition led to expressions for the temperature and chemical potential 
profiles, which turn out to be independent of the nature of the particles, as well as to 
expressions for the electric and heat currents going through the system. Finally, we have 
presented some numerical results, using random matrix theory, supporting the validity of 
Ohm and Fourier laws in our model. 

Let us point out that we could handle the three physical cases, Maxwell-Boltzmann, 
Fermi-Dirac and Bose-Einstein, on the same footing because they share the following two 
features: (i) the specific form (3.4) of / and (ii) the ratio TZ, defined in (3.22), satisfies 
< 7^ < 1. Basically, we expect that any distribution function / satisfying (i)-(ii) will 
lead to the same results. 

For further investigations, it would be interesting to study the general energy-dependent 
S-matrix situation and to compare the corresponding predictions (under RMT) with those 
of the EY-model. One may also analyse the nonlinear transport properties through the full 
counting statistics (ECS) [31,53-59]. More mathematically oriented investigations would 
be to obtain some rigorous results concerning the global composite scattering matrix S and 
consequently to prove the existence of the finite limit, limjv^oo >^{N), which ensures the 
validity of Ohm and Fourier laws in our model. 
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Appendix A: The Coefficients C(n) 

Here we derive some properties of the coefficients C(0), C(l) and C(2). For this, let 
n = 0, 1 or 2. Then 

(i) If xo e (-1, oo), then C^^(n) > 0. Indeed, this follows from 

C^^{0) = e-"° , C^^(l) = (1 + a;o)e-"» , C^^(2) = (2 + 2xo + xl )e-^° . (A.l) 

(ii) If Xq e (— oo, oo), one has C™(n) > 0, and if Xq G (0, oo), one has C^^(n) > 0. 
Indeed, if xq > 0, then clearly C^°(n) > and C^^(n) > 0. For the Fermi-Dirac case. 
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assume that xq < 0, then 

C"°(2) = / -^^dx>0. (A.3) 



(e^ + 1)2 



Now, 



oo 



Observe that ^^(i) > q if = and C™(1) ^ as ^ -oo. Therefore, the 
property C^°(l) > 0, for all < 0, follows from 

AcFD(l)^_ >0. (A.5) 

Appendix B: The Ratio 7^ 

Here, we show that for all xq. 

n = -^= ^ ^-^^ e (0, 1) . (B.l) 
^/Q^ VC(0) • C(2) ^ ' ^ 

Since Qi > and (^2 > 0, one has 7?. > 0, so it only remains to show that TZ <1. 

(i) Let us first consider the Maxwell-Boltzmann case. Recalling the explicit expressions 
(A.l), one obtains {xq e (—1, oo)) 

7^-= ^ + '"° + i <l■ (B.2) 
2 + 2x0 + 

To show that 7?. < 1 in the other two cases, it is convenient to introduce the real functional 
space I/2([a, oo), dx) with scalar product (•, •)„ and norm || • 

/■oo 

(/l,/2)a= / h{x)h{x)dx, ||/||2 = (/,/)^. (B.3) 
J a 

(ii) For the Bose-Einstein case, one has {xq > 0) 

^ ^ T~\\ — IT ' '^^^ ^ I ^ _ .1 ■ (B-4) 

ll5'l|a;o ' 11*^ 9\\xo p -'-I 

Hence, it is sufficient to show that 

llV^^ll^o < llflUo • Ikfi'lUo ■ (B-5) 

We have 

II ^11^0 = {Vxg,Vxg)xo = {9,xg)xo < MUo ■ \\x g\\xo , (B-6) 
where we have used the Cauchy-Schwarz inequality. 
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(iii) Finally, let us consider the Fermi-Dirac case. If xq > 0, then we proceed as in (ii). If 
xo <0, then we observe that 

n< „ '[,^„^"°,, , with g{x) = -4^, (B.7) 
Ikllo-lk^llo |e^ + l| 

and again proceed as in point (ii). 



Appendix C: The Entropy Production 

In this appendix, we derive the expression (3.26) for the entropy production rate ag. How- 
ever, let us first recall how a real positive semi-definite matrix is defined. Let {e^ } C ]R^^+^ 
denote the canonical basis in which (e,, Lej) — Lij, where (•, •) denotes the usual scalar 
product in M^^+^. Then the matrix L is said to be real positive semi-definite if for any 
vector V e M^^^'' one has 

as = {V, LV)^Y1Y1 ^^i^i^o > ■ (C.l) 

i 3 

Let us decompose the matrix elements Lj^ into a symmetric and antisymmetric part: 

(C.2) 

where 

L^ = ^(L,, + L,,) = L^, and L« = 1(L,, - L,,) = -L^ • (C.3) 

Then 

^s = YlJ2^hViV,. (C.4) 

i 3 

This shows that only the symmetric part of the transport matrix L may contribute to the 
entropy production. Observe that the symmetric transport coefficients, L|^, also satisfy the 
relations (3.20) and (3.21). Since we shall only make use of these properties, we see that 
we can assume that the matrix L is symmetric. 

The first idea is to decompose the sum (C.4) into four sums, one associated to each 
block appearing in the matrix L. To simplify the notation, let us consider that the chain 
contains N — 2 quantum dots and then number as 1 and N the left and right reservoirs, re- 
spectively. We shall then write Xi, . . . , Xj^ for the first N components of V and Fi, . . . , Yat 
for its last N components. One may interpret this decomposition as follows: Xi = d/ii/e 
and Yi = 6Ti/T. Recalling that Lij = Lji, one has 

AT TV TV TV TV AT 

= ^ ^ A-,A-, + E i;(4" + ijiv-.v} + E E i« ^is- . (C.5) 

i=l j=l i=l j=l j=l J=l 

Now, using the relations L-j-* = QiLf^^ = hf^ and hf^^ = Q2Lf^ , one obtains 

TV N 

i=l j=l 
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Let Zj = sfW^^j and 7e = Qxj ^fW^- Then 

N N 



1=1 j=i 

Let us first consider the term: 

AT N 

r, = Y,Y.^f^^,x,. (C.8) 

i=\ j=i 

We write 

N N N 

i=l j=l j=l 

Since = 0, one has L^f = - L^f (1 - %) and therefore 

N N N N 

^1 = - E E - '^^.) + E E ^if^i^A^ - ^^j) ■ (C.10) 

i=l j=l 1=1 j=l 

Now comes a trick: 

W TV TV W 

E E - '^^.) = E E - ^^.) ■ (c-ii) 

i=l j=l i=l j=l 

Hence 

AT 

^i = - E ^S? - 2X,X, + X|) . (C.12) 
This relation is well known (see e.g. [35]). Next we consider the second term: 

N N N N N 

r2 = 27e E E ^if^i^j = 2^ E ^u^iZi + 27e E E - ^^j) ■ (c.i3) 

1=1 j=l i=l i=l j=l 

Using the relation = — L^j^ (1 — and the identities 

AT TV TV Af 

E E Lfl^^^Zi{l - 5,,) = E E ^S^^.^.(l - ^ij) ' (C.14) 

i=l j=l i=l j=l 

and 

N N N N 

E E - ^^^) = E E - ^^^) ' (c-15) 
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one can write 

N N 

r2^Y.Yl 4°^^ i^i^j + ^j^i - ^i^i - ^3^3) (1 - ^i.) ■ (c. 16) 

i=l i=\ 

Writing the term in Z^Zj in a similar manner as 7^, one finally obtains 

N 
i < j 

where 

lij — {Xi — Xj)"^ + [Zi — Zj)^ — 2TZCij , Cij — XiZj + XjZi — XiZ^ — XjZj . 
This is the relation (3.26) used in Section 3.3 to prove that ag >0. 

Appendix D: Equilibrium and Non-Equilibrium States 

In this appendix, we prove the equivalent relations presented in Section 3.4. 
To obtain the first one, 

{System is at equilibrium} ■<=^ {li — and Jj = 0, Vi} , 

let us write the no-current condition (i e {L, R, 1, . . . , N}): 

(0) Sf^j , T (1) 



I, . 5;i<«>^ + i«^ = o, (D.l) 

j 

J, = ;^Lg)^ + L<|>^ = 0. (D.2) 



J 



It is convenient to define X,Y e by 



= ^ and = (D.3) 



Then, using the relations L\f = QiL^f and L\f = Q2L\f, and setting Q = Q2/Q1, one 
sees that the conditions (D.1)-(D.2) can be rewritten as follows: 

L(o)(X + QiF) = 0, (D.4) 

L(°)(X + gF) = 0. (D.5) 

Since TZ 1, it follows that Qi 7^ Q. Assume now that all the currents vanish. Then the 
equations (D.4)-(D.5) imply that 

L(°)X = and L^^^y = . (D.6) 
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Since J2j ^j = 0' the equilibrium situation, Xi = Sfi/e and Yi = 6T/T, is a solution. 
Observe now that (— may be interpreted as the generator of an irreducible Markov 
chain and, consequently, the equation L'^'^^V = has a unique (normalised) solution. 
Hence the system is at equilibrium. Finally, the reverse implication is immediate since 
the equilibrium situation satisfies the equations (D.4)-(D.5). 

We next turn to the equivalences in terms of the entropy production: 

{System is at 1 a — and i System is 

equilibrium J * \ out of equilibrium 

Observe that, since cr^ > 0, the two above statements are contrapositive equivalences. 
Therefore, it is sufficient to prove the first. For this we use the expression (C.17) for ag. 

[^] Assume the system is at equilibrium. Then, Xi = ■ ■ ■ = X^ and Z\ — --- — Zn, and 
therefore lij = for all i, j, from which it follows that cr^ = 0. 

[<^] Conversely, assume that as — 0. Then, since (-L^f) < and lij > for all i < j, 
one has lij — for all i < j. Therefore, by (3.27), one deduces that Xi — ■ ■ ■ — Xn and 
Zi — --- — Zj^, i.e. the system is at equilibrium. 




(J, > . 



Appendix E: The Matrix 

Here we prove that the matrix ^ is real positive definite and consequently that l'q^ has 
positive determinant and is therefore invertible. 
Let W em.^,W ^Q, and 

N N 

T ^ {W, Lg^Py) = Y.Y. ^^^^ ■ (E.l) 
i=i j=i 

The idea is to proceed as in the proof of as > 0, but with the identity 

N 

4°' = -iS?-iS-E4°'(i-*«)- (E.2) 

We obtain 

N N 

i,j = l i=l 
i < j 

This shows that L^^ is real positive definite. 
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Now, let Ai, .... Ajv denote the eigenvalues of Lq \ which might be complex since Lq^ 



is not symmetric in general. 

Assume that A is a 
eigenvector. Then 



(1) Assume that A is a real eigenvalue of L^^ and let v E be a normalised associated 

A = X\\v\\ = (v, Xv) = (v, l'^c'^v) > , (E.4) 



where we have used that Lq^ is real positive definite. 

(2) Assume next that A is a complex eigenvalue: A = a + ib, with a.b E M. and b ^ 0. 
Then, since L^q is a real matrix, both A and its complex conjugate A are eigenvalues 
ofL<°'. 

In summary, we have shown that the eigenvalues of L^^ are either positives or come in 
pairs (A, A) with A 7^ 0. Therefore 

det(4°^) = J]Ai >0. (E.5) 



Appendix F: The Profiles 

Here we explain the details in the derivation of the profiles (4.36)-(4.38). From (4.35), one 
sees that the profiles are obtained by solving the generic equation: 

L^X = - ^ D,Xe . (F.l) 

From Appendix E, the matrix can be inverted to give 

X^ - Y^lL^cY'DeXe. (F.2) 

Recalling that (-D^)j = L^jf, one has in components: 



TV 



(R3) 



Now comes a trick: 



N 



k=l 



With this identity, one has 



N N N 

k=i j=i j=i 

N 

= - Y.ii^cr%L^^iXR - Xl) . (R5) 
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Observe that (Lg^ij = L^j = e^C{0)/h Tij and if^ = -e^C{0)/h t^n- Hence, the 
constant e^C(0)//i involved in {[LQ^]~^)ij and in L^*^ simplifies, and we are left with 
profiles independent of the distribution / describing the reservoirs. The temperature and 
chemical potential profiles are obtained by setting Xj = 6Ti/T, with Tj = T + 6Ti, and 
Xi — SfjLi/e, with — SfjLi, respectively. 



Appendix G: An Expression for Ai 

In this appendix, we derive the alternative expression (4.44) for the coefficient Ai. From 
(4.38), one has 

N 

A = J2(^c% tjK ■ (G.l) 

Now, we write 

(rc')., = d^^(-l)'^' det (rc(j, i)) , (G.2) 

so that 

1 ^ 

It remains to expand det (Fc). One has 

N 

det (Fc) = J2(-iy+' det (rc(j, t)) T,, . (G.4) 



Using the relation 



N 



Tji — -FjL - Tj-R - ^ r^fe , (G.5) 



k = l 



one obtains 

N 



det (Fc) = ^(-1)*+^' det (rc(j, i)) [tjj. + ^.-r] - Ri , (G.6) 

where 



= E E (-1)'^' (rc(i, 0) r,fe . (G.7) 

i=i fe = 1 
Finally, using (G.2), one obtains 

N N N 

i?, = det(rc) Yl E(rc')i,r,/c = det(rc) Yl ^^^ = ^- (G.8) 

k = l J=l fe = 1 

fe ^ i fe ^ i 
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Appendix H: The Global Scattering Matrix Sn 



In this appendix, we explain how to build the scattering matrix S of the system made of 
N quantum dots from the N individual scattering matrices S^^\ . . . , S'^^\ where S^''^ is 
the (M + Mfe + M) X (M + M^. + M) unitary scattering matrix associated with the A;-th 
quantum dot. Here, we have set Ml = Mr = M. 

Let us start with the case N = 2. Let S'^^^ and S'^'^^ denote the scattering matrices of the 
first and second quantum dot, respectively, and let S2 denote the corresponding scattering 
matrix for the system made of two quantum dots (see Figure 15). 




Figure 15: Combination of S^^^ and S^'^^ into the composite S-matrix 5*2. 



From Figure 15, one easily sees that these matrices are related as follows: 





5(2) 





S, 



(xA 






Xl 




yi 








\xrJ 







Introducing 
and 



(1) c(i) 



^11 

'^21 



^12 



.'(1) 



Xhl 



c(i) +(1) 
^33 ' 



0(1) c 
^31 ^' 



(1) 
32 



XL 
Xl 



yLi 



one can rewrite the first relation in (H.l) as follows: 

^r(i) t'(i)\ fxLi 



yhi 



Similarly, setting 

r-{2) _ C:(2) 



' q(2) 
•^22 
c(2) 
•^32 



0(2)' 
'-'23 
C(2) 
"^33 



(2) 



•^21 

c(2) 

'^31 



,t' 



'(1) 



' q(1)' 
"^13 

•^23 



/(2) 



0(2) 
^12 



(H.1) 



(H.2) 



(H.3) 



(H.4) 



(H.5) 
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and 



X2R 



X2 
Xr 



one has for the second relation in (H.l): 

■^(2) f(2) 
t{2) ^'(2) 



Z 
•«2R 



1/2 
2/R 



Z 
2/2R 



(H.6) 



(H.7) 



Finally, we write the third relation in (H. 1) as follows: 



r2 ^2^ (^hi 

t2 rU \X2R 



Vli 

y2R 



(H.8) 



Solving the equations (H.4) and (H.7) for and y2R, and recalling (H.8), one obtains 



r2 = r«+t'(i)[l_r(V«]-V(2)t(i) , 
= /(2)^^(2)[^_^/(i)^(2)]-y(i)^/(2) ^ 

4 = 



^(2)/(l)l-l^'(2) 



(H.9) 
(H.IO) 
(H.ll) 
(H.12) 



Let us now generalise these relations to the > 2 quantum dots situation. For this we 
denote by S^^\ . . . , S^^^ the N individual scattering matrices and by Sn the composite 
S-matrix associated to N quantum dots. We write 



S 



N 









Xi 




yi 






yN 


\Xr/ 




WJ 



s 



N 



(H.13) 



Then 



j.^ ^ ^(12...7V-l) ^ ^/(12...iV-l)j^ _ ^(iV)^/(12...iV-l)j-l^(iV)^(12...iV-l) 

j.'^ ^ ^/(JV) _^ ^(7V) _ ^/(12...iV-l)^(7V)|-l^/(12...iV-l)^/(JV) ^ 

tj, = i(iV)[i_/(12...iV-l)^(JV)]-l^(12...7V-l)^ 

4 = i'(12...iV-l)[i_^ 



,(JV)^/(12...7V-l)i-l^/(JV) 



(H.14) 
(H.15) 
(H.16) 
(H.17) 



where 

rim 

and 



^22 ^23 1 JN) 
^32 '-'33 




AN) 
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(H.18) 



,(i2...Ar-i) 



'(-^iV-l) 



11 



{Sn-i) 



IN 



,(5'iv-i)ivi ■ ■ ■ (Sn-i) 



, r'(^'-^-^) = (5^-i)iv+i iv+i , (H.19) 
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^{12. ..N 1) _ (J^Sn-i)n+1 1 • • • {Sn-i)n+1 n) , t' 



/{12...N-1) 



1 N+l 



N-1)N 7V+1/ 



(H.20) 

Therefore, by working recursively, one can build the global scattering matrix S = Sn from 
the local scattering matrices S^-^\ . . . , S^^K Note that, by construction, if S^^\ . . . , S^^'> 
are unitary, then S will also be unitary. 

Remark. If one considers that the terminals have only one channel: Mj = 1 for all i, then 
r^^^ and r'^^^-^^ become complex numbers and we can introduce the following complex 
number: 

c^= [l-rWr'(i2-^-i)]-^ (H.21) 
Then one can rewrite the equations (H.14)-(H.17) as follows: 

= ^(12...7V-l) ^ ^(iV) ^/(12...iV-l)^(12...iV-l) ^ (H_22) 

r'^ = r'(^) + c^r'(i2...7v-i)^(iv)^/(iv)^ (H_23) 

= CArt(^H(i2-^-i), (H.24) 
4 = c^t'^^^-^-'h'^''K (H.25) 

These are the equations that are used in the numerical simulations in Section 5. 
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